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1. Introduction 



A scheme for investigating amplification and discrimination of input photons is presented 
based on the interaction of a large number of identical two-level atoms and low intensity 
single mode radiation. This analysis provides a deeper understanding of the cumulative 
response of millions of atoms to few input photons within the basic model [1] of atom- 
field interaction. 

Historically, collective behavior of iV two level atoms with single mode radiation, 
investigated by Dicke [1] in 1954, has led to vast array of interesting physical phenomena 
in quantum optics and recently, in artificial condensed matter systems [2]. The present 
investigation is an entirely novel application of the Dicke model. The fact that the sum 
of the number of photons and the number of atoms in the excited states is a constant 
under time evolution has a natural consequence of rearranging the atom-photon numbers 
resulting in photon amplification as well as discrimination of small number of photons. A 
projective measurement on temporally evolving combined multiatom-radiation system 
is shown to lead to both these features. 



2. Physical Model 

The Hamiltonian characterizing the interaction of N two level atoms with a single mode 
radiation is given by (H = 1) 

H = uj(o)a) + uj S z H — -= (S+ a + S- afi), (1) 

v N 



where g denotes the atom-photon coupling parameter; u , the atomic splitting and u 
the filed frequency; a* (a) are creation (annihilation) operators of the field satisfying 
the bosonic commutation relations [a, a'] = 1. The collective (pseudo) spin operators 
of the two-state atoms 

^ N ^ N 

S± = S x ±iS y = - a a±, S z = - ^2 a az (2) 

^ a=l ^ a=\ 

obey the commutation relations 

[S + ,S^] = 2S Z , [S Z ,S±)=±S±. (3) 
It is well-known [3] that the excitation operator 

N cx = a)a + S z + — (4) 

remains constant as the system evolves under the Hamiltonian (TjQ). The degenerate 
atom-photon states |S'= ; j,M = n e — ® |w) (— y < M < or equivalently 
n e = 0, 1, 2, . . . , N, and n = 0,1,2,...), with a fixed eigenvalue [j n + n e of (BJ span 
the space of the Hamiltonian. The symmetric atomic Dicke states IS* = y, M = n e — y) 



% The Fock states \n), n = 0,1,2,... denote eigenstates of the photon number operator a^a with 
eigenvalue n and the Dicke states \S = =j, M), — y < M < ^ are joint eigenstates of S 2 = S 2 + Sy + S 2 
and S z with eigenvalues y(t + ■'■)' ^ res P ec tively. 
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(except for n e = and N) are well-known for their entanglement properties [U EJ E|. 
These are essential in the ensuing discussion The collective ground state of the atomic 
system is denoted by |0) = 1^, — y) throughout this sequel. 

A given initial state of the combined atom-radiation system, Par(O) evolves to 

Par® = e- iHt PAR(0)e iHt . (5) 

Consequently, a projection operator L~Uo <8> Ir, with ITao = |0)(0|, and Ir the unit 
operator in the radiation space, gives us the conditional density matrix of the collective 
atom-photon system ((where all the atoms are projected to their collective ground 
state)), subjected to the constraint (TjJ. This projective measurement has the implication 
that maximum number of photons are emitted consequently. The collective atomic 
ground state serves as a detection device and the projection operation at subsequent 
suitable intervals of time corresponds to efficient measurement of the radiation state 
containing maximum number of photons allowed. It may be observed that the constant 
of motion leads to rearrangement of the number of atoms in the excited state with the 
number of photons, leading to photon amplification. This is evident if n e > n number of 
atoms are initially in excited state with n input photons, in which case (N ex ) = n+n e and 
a subesquent measurement projects the atoms to their collective ground state leading 
to (a^a)(t) = n + n e . More generally, (aJa)(t) = (iV ex ) — (S z )(t) — y. A judicious 
choice of the initial state and the time of projection measurement is the basis of photon 
amplification in this model. When a pure collective low lying excited state of atoms are 
considered, the method outlined above results in photon number discrimination as well. 



3. Photon amplification with different initial states 

3.1. Pure atom-photon state \n e ;n) 

Let us consider an initial atom-photon state 
N N 

\S = — ,M = n e - —) ®\n) = |n e ;n), (6) 

with n e « N, and the number of atoms iV sufficiently large so that the Holstein- 
Primakoff mapping [7j 

N 



S + = tf^N- tfb, S_ = yjN - tfb b, S z = tfb-— (7) 

in terms of bosonic operators b, satisfying 

[6,6t]=l, (8) 

reduces the Hamiltonian ([T]) into a two mode bosonic interaction structure [2] (up to 
order 0(1/ N)): 

N 

H ~luN cx -lu — + g(b j a + ba j ). (9) 



§ Dicke states with a specific permutation symmetry are chosen here because of the greatest simplicity 
offered by them in their theoretical analysis and also due to their current experimental relevance [H [S] . 
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Temporal evolution in this approximation thus assumes the form, 

U(t) = e~ iHt = e -*Ti? N --%*O e -iT<bta+*afi) m (10) 

(Here we have denoted gt = r.) It may be noted that U(r) acts as a passive unitary 
beam-splitter on the two mode bosonic Fock states |n e ;n). Thus, initial states of the 
form get confined within the space spanned by the (n e + n + l) basis states \n' e ; n'), 
with n' e , n' — 0, 1, 2 . . . such that n' e + n' = n e + n under time evolution. It is interesting 
to note that at scaled time r = n/2 the unitary operator U(t) swaps [8] the atom-photon 
numbers i.e., \n e ;n) — > \n;n e ) with unit probability. 
In general, we have, 

U(r)\n e] n) = e- ir[ ^^ + ^ ] £ \n' e ,n')(n' e ,n'\e- iT ^ a+ba ^\n e ,n) 

■ r (n + n e )uj , Nun n ri e +n 

= e~^— + ^ £ \n' e ,n'}d n ,*,(2T) 



where [9] 



^ (i + m )Ki - m )-U + m 'Y-U - m 'V- , 



fc— m+m' 



rf m', m ( 2r ) E ^ + m + ^ !A .,^ _ k _ m ,y^ k _ m + m y 

x (cosr) 2 ^ 2fe+m - m ' (sinr) 2fe ~ m+m ' = ^ m ,(-2r) (12) 

with the sum over k taken such that none of the arguments of the factorials in the 
denominator are negative. 

A measurement Il^o <8> Ir at an instant t projects all the atoms to ground state: 

n e -\-n 

n A o ®IrU (t) |n e ; n) = J] |0, n') d f n+ne) „ e _„ (2r) £ n /, ne+n 

2 ' 2 

3J7 r[(ne-n+n')]/4 -it [ <■ + e > 



e 



-it{nu + ^ ) g-ineC^-f) (.^e ^v[n e ,n,T) |0; n + n e ).(13) 

Thus, the probability of finding the atoms in ground state, which in turn corresponds 
to that of photon amplification n — > n + n e , is given by 

V(n e , n, r) = ^ n + ne J cos 2 "(r) sin 2 " e (r). (14) 

In Figs. [1] (a), (b), (c), we have displayed the probability of maximum photon emission 
for three choices of n e = 1, 10, and 25, each with different values for the number of 
photons n = 0,1,5,10, as a function of scaled time r. For n = i.e., dark photon 
input state, the probability has a single peak [jj] in the range < r < ix with maximum 
value V(n e ,0,T = vr/2) = 1, which corresponds to swapping of atom-photon numbers 
|n e ; 0} — > |0;n e ). However, when n ^ the projective measurement n^ ® Ir leads 
to \n e ;n) — > |0;n + n e ) corresponding to maximum photon emission; this is not a 

|| The probability "P(n e ,n, t) exhibits repetitive structure with a period 7r, as is evident from l|14[). 
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Figure 1. Probability V(n e , n, r) of finding the atoms in ground state as a function of 
scaled time r for different values of initial photon numbers n. Maxima in probabilities 
occur at scaled times t„ — arccos , / — £ — • 



complete swap action and occurs with probability < V(n e ,n ^ 0, r) < 1. Clearly, 
the probability V(n e ,n ^ 0, r) vanishes at scaled time r = n/2, in order to give 
way to complete swap action \n e ;n) — > \n e ;n). The probability profile for maximum 
photon emission n — > n + n e therefore reveals two peaks around r = ir/2 for n / 0. 
The maxima of probabilities corresponding to different input photon numbers are well 
separated and they appear periodically at r n = arccos y^j^-, allowing for photon 
number discrimination. 

We define the time of perception T p (n e ,n), for a given n e and n, as the time at 
which the probability V(n e , n, r) attains its maximum. This determines the efficient 
detection of the photons. An examination of these figures reveal: (i) for a given n e , 
the time of perception r p (n,n e ) reduces as more and more photons are detected - with 
lesser and lesser efficiency. Moreover the widths in the probabilities V(n e ,n,r) reduce 
correspondingly, (ii) As n e increases, there is a threshold time for the detection of 
photons, before which the probabilities are zero. In particular, for n e = 1, there is an 
instant response to photons (for all n) as seen in Fig. [T](a), whereas for higher values 
of n e there is a delay in such a response (See Fig. [D (b), and (c)). (iii) For a given 
photon number n the profile of probability as a function of time sharpens as n e increases; 
however maximum value of the probability drops with this, as is evident from Figs. HJ 
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3.2. Pure state of atoms with low intensity coherent radiation 

The above discussion was confined to pure photon number states. We now show that 
an enhanced photon amplification behavior is realized, when initially a low intensity 
coherent state of radiation 

oo n 

\a) = e'\ a \ 2 ' 2 f] ^=\n), \a\ 2 <\ (15) 

is considered as input. The probability of finding the atoms in ground state is obtained 
by following the procedures given in (JTTj) and (|T3l) : 

V(n e ,a,r) = e"H 2 £ p( ne>n , r )L_L_. (16) 




Figure 2. Probability V(n e ,a,T) of finding the atoms in ground state as a function 
of scaled time r for different values of initial intensity \a\ 2 and for different choices of 
n e , the initial number of atoms in excited state. 

The probability profile f|T6|) - with coherent radiation as input - is a series involving 
individual atom-photon Fock state probabilities V(n e ,n,r), dark photon probability 
V(n e , 0, r) being the leading term for low value of intensity |a| 2 . Thus, the probability 
response to low intensity coherent light has a similar structure as that of dark photon 
probability V(n e , 0, r) (see Fig. [[]) with small contribution from higher order terms 

|a| 2n , n = 1, 2, It may be seen from Fig. [2] that the probability V(n e , a, r) is found 

to be nearly zero in the begining and raises to a maximum at r = tt/2 (peak value 
being V(n e ,a,n/2) ~ e" |a|2 ). An increase in the initial intensity of coherent radiation 
has the effect of reducing the probability of finding the atoms in their ground state. 
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Figure 3. Probability of finding the atoms in ground state under time evolution 
governed by the Hamiltonian 0, when initial atomic state is chosen as (i) a mixed 
state - 1 +1 X)m=o l m )( TO l; ( u ) P ure state \n e ), and the radiation in a coherent state, 
| a), as a function of scaled time r, for two different values of inital intensity of radiation 
and n e = 25. 



Amplification of the intensity of radiation, j^gj = i °ja , (after performing projective 
measurement Uao on the temporally evolving state U{r)\n e ] a)) approaches the value 
l^p as r — ► | i.e., when the probability of finding the atoms in ground state is maximum. 

3.3. Mixed state of atoms with low intensity coherent radiation 

So far, we have considered pure atomic states. It is of interest to investigate how the 
effects found above may get modified when mixed atomic states are employed. For 
simplicity, we choose here a chaotic mixture (n e + 1) of low lying collective excited 
states of atoms, 

i n e 

Patom = — r V] \m)(m\, n e «N, (17) 

U e + 1 ^0 

along with low intensity coherent radiation. The probability of finding the atoms in 
ground state under time evolution in this model is readily found to be, 

e -M 2 n <= 00 \ a \2n 

"P(Patom,a,r) = — ^ X — r V(m,n,T). (18) 

The probability V(p at om, a, t) as well as V(n e , a, r) (see Eq. (fT6l) ) associated with a pure 
atomic state |n e ) are plotted, as a function of scaled time r, in Fig. [3] (a) and (b) for 
two different choices of the initial intensity of radiation for a fixed n e = 25. We observe 
that while the maximum probability remains the same for a given initial intensity of 
radiation \a\ 2 , the mixed states p a tom result in a wider spread around the maximum 
value (at r = 7r/2) compared to their pure state counterparts. Also, the probability 
^(Patom, ot, t) builds above the background value ^rj. Moreover, as the intensity of the 
radiation \a\ 2 increases, the maximum probability drops down. It may be noted that 
the photon amplification factor approaches the value 2 ^ 2 as r — ► n/2 for mixed 
atomic states p at0 m m contrast to its corresponding value r^fe in the case of pure states 
\n e ). 
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4. Summary 

We have presented a scheme for photon amplification and discrimination based on 
the Dicke model interaction between single mode radiation and iV-atom system. This 
involves Holstein-Primakopf large N approximation [Tj, where the number n e of atoms in 
excited states is assumed to be much smaller than the total number N of atoms. Further, 
this approximation enforces that the number n of input photons is small, compared 
to N. This leads to a beam splitter feature for time evolution under the interaction 
Hamiltonian (Q, thus ensuring entanglement between the atoms in collective excited 
states and the photons [TU] . Emission of maximum number of photons corresponds 
to a completely uncorrelated atom-radiation system, with all the atoms in the ground 
state and this is acheived via a projective measurement at a suitable interval of time. 
Concepts like threshold time and time of perception for exposure to low intensity light 
emerge as characteristic features of our investigation. Discrimination of small number 
of photons is also realized in this scheme. This provides motivation for a quantum 
mechanical analysis of the collective response of millions of rods in the eye, which act 
as nearly perfect photon detectors, initiating the associated process of vision during 

night PHE21H3]. 
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